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Abstract. Let fi,...,f a G K[xi, . . . ,x m ] be a system of polynomials gener- 
ating a zero-dimensional ideal X, where K is an arbitrary algebraically closed 
field. We study the computation of "matrices of traces" for the factor algebra 
A := K[zi, . . . , x m ]/X, i.e. matrices with entries which are trace functions of 
the roots of I. Such matrices of traces in turn allow us to compute a system 
of multiplication matrices {M Xi \i = 1, . . . , m} of the radical 

Wc first propose a method using Macaulay type resultant matrices of /i , . . . , f s 
and a polynomial J to compute moment matrices, and in particular matrices 
of traces for A. Here J is a polynomial generalizing the Jacobian. We prove 
bounds on the degrees needed for the Macaulay matrix in the case when I has 
finitely many projective roots in PJ^. We also extend previous results which 
work only for the case where A is Gorenstein to the non-Gorcnstcin case. 

The second proposed method uses Bezoutian matrices to compute matrices 
of traces of A. Here we need the assumption that s = m and f\ , . . . , f m define 
an affine complete intersection. This second method also works if we have 
higher dimensional components at infinity. A new explicit description of the 
generators of \/X are given in terms of Bczoutians. 

Keyword: matrix of traces; radical of an ideal; 

1. Introduction 

This paper is a continuation of our previous investigation in [251 126] to compute 
the approximate radical of a zero dimensional ideal which has zero clusters. It turns 
out that the computationally most expensive part of the method in [25l [26] is the 
computation of the matrix of traces. We address this problem in the present paper. 
Some of the results of this paper also appeared earlier in [23], however here we 
present generalized versions of those results and also add new results, as described 
below. 

The computation of the radical of a zero dimensional ideal is a very important 
problem in computer algebra since a lot of the algorithms for solving polynomial 
systems with finitely many solutions need to start with a radical ideal. This is also 
the case in many numerical approaches, where Newton-like methods are used. From 
a symbolic-numeric perspective, when we are dealing with approximate polynomi- 
als, the zero-clusters create great numerical instability, which can be eliminated by 
computing the approximate radical. 

The theoretical basis of the symbolic- numeric algorithm presented in [25]. [26] was 
Dickson's lemma |15) . which, in the exact case, reduces the problem of computing 
the radical of a zero dimensional ideal to the computation of the nullspace of the 
so called matrices of traces (sec Definition I3.17j) : in [25J [26] we studied numerical 



This research was partly supported by the Marie-Curie Initial Training Network SAGA. 

1 



2 I. JANOVITZ-FREIREICH, B. MOURRAIN, L. RONYAI, AND A. SZANTO 

properties of the matrix of traces when the roots are not multiple roots, but form 
small clusters. Among other things we showed that the direct computation of 
the matrix of traces (without the computation of the multiplication matrices) is 
preferable since the matrix of traces is continuous with respect to root perturbations 
around multiplicities while multiplication matrices are generally not. 

In the present paper, first we give a simple algorithm using only Macaulay type 
resultant matrices and elementary linear algebra to compute matrices of traces of 
zero dimensional ideals which have finitely many projective roots. We also extend 
the method presented in [24] to handle systems which might have roots at infinity 
or for which the quotient algebra is non-Gorenstein. 

In the second part of the paper, we investigate how to compute matrices of traces 
using Bezoutians in the affinc complete intersection case. Our approach in that case 
is based on [40l [39] . 

For the method using Macaulay matrices we need the following assumptions: let 
f = [/i, . . . , f s ] be a system of polynomials of degrees d\ > ■ ■ ■ > d s in K[x], with 
x = [x\ , . . . , x m ] , generating an ideal X in K[x] , where IK is an arbitrary algebraically 
closed field. We assume that the algebra A := K[x]/X is finite dimensional over K 
and that we have bounds 8 > and < k < 8 such that a basis B = [6i, . . . , 6jy] of 
A can be obtained by taking a linear basis of the space 

W(/i,..,/ s ) s nK[x] t 

where K[x]fc is the set of polynomials of degree at most k and (/i, . . . , f s )g = 
{J2i=i1ifi '■ degg.; < 8 — di}. We can assume that the basis B consists of 
monomials of degrees at most A: by a slight abuse of notation. In our earlier work 
[24j . we gave bounds for k and 8 in the case where there were no roots at infinity 
using a result of Lazard [37] (see Theorem l3.4p . Here we extend those results to the 
case where X has finitely many projective common roots in PJ^ (see Theorem 13.51) . 
Furthermore, we now extend the method presented in |24| . which only addressed the 
case where A is Gorenstein over K (see Definition 13. 1| . to handle non-Gorenstein 
algebras. 

The main ingredient of our first method is a Macaulay type resultant matrix 
MacA(f), which is defined from the transpose matrix of the degree A Sylvester 
map (<?i, . . . ,g s ) i— > X);=i fidi G IK[ x ]a for A < 28 + 1 using simple linear alge- 
bra (see Definition 13. 8|) . Using our results, we can compute a basis B of A using 
MacA(f)- We also prove that a random element y of the nullspace of MacA(f) pro- 
vides an N x N moment matrix 9JIb (y) which has the maximal possible rank with 
high probability (similarly as in [35] ) . Note that in the Gorenstein case the moment 
matrix 97ls (y) is non-singular. This will no longer be true in the non-Gorenstcin 
case. This moment matrix allows us to compute the other main ingredient of our 
algorithm, a polynomial J of degree at most 8, such that J is the generalization of 
the Jacobian of f±, . . . f s in the case when s = m. The main result of the paper 
now can be formulated as follows: 

Theorem Let B = [bi, . . . , fejv] be a basis of A with deg(6i) < k. With J as above, 
let Sylg(J) be the transpose matrix of the map YliLi l— * J • Y^iLi G K[x]a 
for Ci £ K. Then 

[T 1 r(6 J :6 J )]5 =1 =Syl s (J)-X, 
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where X is the unique extension of the matrix 971s (y) such that MacA(f) • X = 0. 

Once we compute the matrix of traces R := [Tr(bibj)]f^ =1 and the matrices 

R Xk := [Tr(xkbibj)]fj =1 = Syl B (xkJ) ■ X for k = 1, . . . ,m, we can use the results 
of [35[ [55] to compute a system of multiplication matrices for the (approximate) 
radical of X as follows: if R is a (numerical) maximal non-singular submatrix of R 
and R Xk is the submatrix of R Xk with the same row and column indices as in i?, 
then the solution M Xk of the linear matrix equation 

RM Xk = R Xk 

is an (approximate) multiplication matrix of Xk for the (approximate) radical of 
X. See [26j for the definition of (approximate) multiplication matrices. Note that 
a generating set for the radical VX can be obtained directly from the definition of 
multiplication matrices, in particular, it corresponds to the rows of the matrices 
M Xl , . . . , M Xm . 

We also point out that in the s = m case these multiplication matrices M Xk 
of Vx can be obtained even more simply using the nullspace of MacA(f) and the 
Jacobian J of f , without computing the matrices of traces. 

In the last section we investigate the use of Bczoutians to compute matrices of 
traces of systems fi, ■ ■ ■ ,f m which form an affinc complete intersection. In this 
particular setting, our method allows systems that may have higher dimensional 
projective components. 

In the univariate case it is proved in |40j that the Bezoutian matrix of a univariate 
polynomial / and its derivative /' is a matrix of traces with respect to the Horner 
basis of / (see subsection 14. lj) . Therefore, applying the method to compute the 
approximate or exact radical from the matrix of traces provided by the Bezoutian 
will give us an approximate or exact square-free factorization of /. The question 
that naturally arises is how this method relates to computing the square-free factor 
as gcd (f ji\ ■ We show here that the two algorithms are computationally equivalent. 

The generalization to the multivariate case is not quite as straightforward. The 
goal would be to express the Bezout matrix of f\ , . . . , f m and their Jacobian J as 
a matrix of traces with respect to some basis, generalizing the univariate case (see 
the definition of the Bezout matrix - sometimes also referred as the Dixon matrix 
- in Definition ^. 5p . Unfortunately, the Bezout matrix cannot directly be expressed 
as a matrix of traces. However, in [40] it is shown that a reduced version of the 
Bezout matrix of fi , . . . , f m , and J is equal to the matrix of traces of f\ , . . . , f m 
with respect to the so called canonical basis, obtained from the reduced Bezout 
matrix of fi , . . . , f m , and 1 . The required reduction of the Bezout matrix involves 
reducing polynomials modulo X. 

Now the question is how to find the reduced version of the Bezout matrix without 
further information on the structure of the quotient algebra C[ari, . . . , x m ]/I, e.g. 
without Grobncr Bases or multiplication matrices. First we show that we can obtain 
a set of generating polynomials for the radical y/l from the non-reduced Bezoutian 
matrices (see Theorem 14. 9p . Secondly, we give an algorithm which computes a 
system of multiplication matrices M Xl , . . . , M Xm for \fl. This algorithm adapts the 
results of [39] to find the required reduced Bezout matrices using only elements in 
\fl which were obtained from non-reduced Bezout matrices. 
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2. Related Work 

The motivation for this work was the papers |35[ 136] where they use moment 
matrices to compute the radical of real and complex ideals. They present two 
versions of the method for the complex case: first, in [36] they double up the 
machinery for the real case to obtain the radical of the complex ideal. However, in 
[35] they significantly simplify their method and show how to use moment matrices 
of maximal rank to compute the multiplication matrices of an ideal between X 
and its radical y/T. In particular, in the Gorenstein case they can compute the 
multiplication matrices of X. In fact, in |35| they cite our previous work [25] to 
compute the multiplication matrices of \/X from the multiplication matrices of X, 
but the method proposed in the present paper is much simpler and more direct. 

Note that one can also obtain the multiplication matrices of X with respect to 
the basis B — [bi, . . . , 6 at] by simply eliminating the terms not in B from Xkh using 
Mac,5+i(f). The advantage of computing multiplication matrices of the radical y/T 
is that it returns matrices which are always simultaneously diagonalizablc, and 
possibly smaller than the multiplication matrices of X, hence easier to work with. 
Moreover, if B contains the monomials 1, x\, . . . , x m , one eigenvector computation 
yields directly the coordinates of the roots. 

Computation of the radical of zero dimensional complex ideals is very well studied 
in the literature: methods most related to ours include [30] [3] where matrices 
of traces are used in order to find generators of the radical, and the matrices of 
traces are computed using Grobner Bases; also, in [TJ they use the traces to give a 
bound for the degree of the generators of the radical and use linear solving methods 
from there; in [3TJ they describe the computation of the radical using symmetric 
functions which are related to traces. One of the most commonly quoted method 
to compute radicals is to compute the projections X n K[a;,] for each i = 1, . . . , m 
and then use univariate squarefree factorization (see for example [TU [30] [Til HI] 
). The advantage of the latter is that it can be generalized for higher dimensional 
ideals (see for example [29]). We note here that an advantage of the method using 
matrices of traces is that it behaves stably under perturbation of the roots of the 
input system, as was proved in [26] . Other methods to compute the radical of zero 
dimensional ideals include [28] [TU [32] [33] [34] [45] . Applications of computing the 
radical include [23], where they show how to compute the multiplicity structure of 
the roots of X once the radical is computed. 

Methods for computing the matrix of traces directly from the generating poly- 
nomials of X, without using multiplication matrices, include [TU [S] where they use 
Newton Sums, [8] [9] [10] where they use residues and [13] using resultants. Besides 
computing the radical of an ideal, matrices of traces have numerous applications 
mainly in real algebraic geometry [4] [41] [5] , or in [42] where trace matrices are 
applied to find separating linear forms deterministically. 

3. Ideals with Finitely Many Projective Roots 

3.1. The Gorenstein Case. Some of the results of this subsection appeared in 
[24] . We included them here for completeness. 

Let f = [fi , . . . , f s ] be a system of polynomials of degrees d\ > ■ ■ ■ > d s in 
K[x], where x = [xi, ... ,x m ] and K is an arbitrary algebraically closed field. Let 
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X be the ideal generated by f±, . . . , f s in K[x] and define A := K[x]/X. We assume 
throughout the paper that A is a finite dimensional vector space over K and let A* 
denote the dual space of A. 

Let us first recall the definition of a Gorenstein algebra (c.f. [311 [43j HZ1 [35]). 
Note that these algebras are also referred to as Frobenius in the literature, see for 
example [2]- 

Definition 3.1. A finite dimensional ¥L-algebra A is Gorenstein (over KJ if there 
exists a nondegenerate 'K-bilinear form B(x,y) on A such that 

B(ab, c) = B(a, be) for every a, 6, c £ A. 

Note that this is equivalent to the fact that A and A* are isomorphic as A 
modules. It is also equivalent to the existence of a K-linear function A : A — > K 
such that the bilinear form B(a, b) := A(ab) is nondegenerate on A. 

Assumption 3.2. Throughout this subsection we assume that A is Gorenstein. 
Furthermore, we also assume that we have a bound S > and < k < S such that 



(1) dim K K[x] fe /(/i, . . . , f a ) s n K[x] fc = dim K K[x] fc /(/i, . . . , f s ) d n K[x] fc 
for all d> S. Here K[x]fe := {p e K[x] : deg(p) < k} and 

(2) (fu f s )d ■■= fe /** : d °s(*) < d - *| • 

Theorem 3.3. Assume that 8 and k satisfy the condition {!]]. Then 
dim K (A) = dim K K[x] fc /(/ 1 , ...,/»>* n K[x] fc . 

Proof. Assume that S and k satisfy the condition ([1]) and let B := [i>i , . . . , £>jv] 
be a basis for K[x]j,/(/ l5 . . . , / s )^ n K[x]fc. Taking prc-images, we can assume 
that bjv arc polynomials in K[x]fc. We claim that B is a basis for A = 

K [ x ]/(/i> •••,/«)• Since . . . , f s )s Q (f ± , f s ) d if 6 < d, B is clearly a gener- 
ator set for A. On the other hand, assume that B is not linearly independent in 
A, i.e. there exist c\, . . . , cjv S K such that Xa=i ^ s m (/i) ■ ■ ■ i /«)• Then there 
exists d > 8 such that c,^ G . . . , / s )d- But J^i=i c s^i i s a ^ so m IK [x] ^ , so 

B is linearly dependent in K[x]fc/(/i, . . . , f s )d n K[x]fe, which contradicts condition 
©. □ 

We have the following theorems giving bounds for 5 in the case when f has 
finitely many projective roots. First we assume that f has no roots at infinity. 

Theorem 3.4. Let f = [fi,---,f s ] be a system of polynomials of degrees d\ > 
■ ■ ■ > d s in K[x] . Assume that the corresponding system of homogenous polynomials 
fi, ... ,fg has finitely many projective common roots in P^ 1 . Assume further that 
fx , . . . , f s have no common roots at infinity. Then: 

(1) If s = m then for S — k := X^Li (di — 1) condition (QJ) is satisfied. Further- 
more, in this case A is always Gorenstein. 

(2) If s > m then for 5 = k := X^i 1 ~ 111 condition (QJ) is satisfied. 

Proof. For the first assertion let f h = [f^, . . . , /^] be the homogenization of f using 
a new variable x m +\. Using our assumption that f h has finitely many roots in PJ^ 1 
and s = m, one can see that (i h ) is a regular sequence in R := K[xi, . . . , x m ,x m +i]. 
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Define the graded ring B := R/(i h ). Following the approach and notation in [44], 
we can now calculate the Hilbert series of B, defined by H(B,X) = J2d T~(-B(d)\ d , 
where Hb is the Hilbert function of B. We have 

H(R > X)= H ™ 



(1_ \di) ...(i_ A d ™)' 
and using the simple fact that 

H(R, A) 



(1 - A) m +! 
we obtain that 

(l + A + .-. + A^-^-.-ll + A+.-. + A^- 1 ) 

= 

= 5(A)(1 + A +...), 

where 

5 (A) = (1 + A + • • • + A^- 1 ) • • • (1 + A + • • • + A 11 '"- 1 ). 
This implies that the Hilbert function 

H B (8) =H b (6 + 1) = H b (5 + 2) = ... 

Note that dehomogenization induces a linear isomorphism Bd — ► K[x]d/ • . • , f m )d, 
where Bd stands for the degree d homogeneous part of B. From this, using that 
there are no common roots at infinity, wc infer that for d > 5 dim^ K[x]^/ . . . , f m )d 
dimjc^l = N, which implies (Q]). 

Note that the common value N = Hb(S) is the sum of the coefficients of g, which 

is 

m 
i=l 

To prove that A is Gorenstein, we cite [17j Proposition 8.25, p. 221] where it 
is proved that if fx, . . . , f m is an affine complete intersection then the Bezoutian 
B\.f 1 ....j m defines an isomorphism between A* and A. 

To prove the second assertion we note that [371 Theorem 3.3] implies that 

dim K B s = dim K B s+ i = 

From here we obtain (p} as in the Case 1. □ 

The following theorem generalizes the previous result for systems which may 
have roots at infinity. 

Theorem 3.5. Let f = • • • , /s] be a system of polynomials of degrees d\ > 
■ ■ ■ > d s in K[x]. Assume that the corresponding system of homogenous polyno- 
mials fi,...,fs has finitely many projective common roots in PJj? and let J := 
(/■L , . . . , /g ) be the ideal they generate in R := K[xi, . . . , x n , x n +i]. Then: 

(1) If s — m then for k := (dj — 1) and S := k+l condition (Qp is satisfied. 

(2) If s > m then for k := ^"^j 1 dj — m and 5 := k+1 condition (QJ) is satisfied. 

Proof. Assume that f has N affine roots and N' roots at infinity, counted with 
multiplicity. In this proof only, for a homogeneous ideal J C R, J t denotes the 
elements of J of degree equal to t, abusing the notation. By the proof of Theorem 
and [37], we have that for k defined above (in both cases) and for all d > 0, 
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Rk+d/Jk+d (resp. Rk+d/{J + (x m +i))k+d) is of dimension JV + JV' (resp. N'). 
Consider the exact sequence 

— > ( J : x m +i)k/Jk — > Rk/Jk — Rk+l/Jk+l — > Rk+i/(J + x m +i)k+i — > 

where M m+ i is the multiplication by a; m+1 . Using the relation on the dimensions 
of the vector spaces of this exact sequence, we deduce that 

dim K ((J : awiWJfc) = (JV + JV') - (JV + JV') + JV' = JV'. 

Thus we can choose a basis \b\, . . . , frjv+w 

] of i? fe / Jfc such that G (J : 

2Wfi)k- Moreover, we can even assume that b^, +1 , . . . , b^ +N , G (z m +i)k> since 

span(&J, . . . ,b% +N ,) = span(foi, . . . ,b%,) + {x m +i)k + Jk- 

If x^ n+1 G {J : x rn +i)k then xj^+j G J, all the roots are at infinity, JV = 0, and 
{•J '■ x m +i)k = Rk modulo Jk which shows that x m +iRk C J7fc+i- After dehomo- 
geneization all polynomials of degree < k are in (/i, . . . , f s )k+i- So condition {l} is 
satisfied for 5 = k + 1. 

Suppose now that G" (J7" : x m +i)k, so that we can take b^, +1 = x^ l+1 . 

As x m +i(J : x m +i) G J , we deduce that x m +i fe'j; = • • • = a; m +i 6^/ = modulo 
J7a:+i and that 

dim K (span (x m +i &i, •• - ^m+i &/v+jv) /Jk+i) < N. 

As we have 

= span {b\, . . . , b h N+N ,) + J h 
and x m+ \ b\ G Jk+i for 1 < i < JV', we deduce that 

Xm+l Rk = Span (x m -)-l6^r/ + 1 , . . . , .T m +l6^r + jY/ 

After dehomogencization, we obtain a family B = [£>jy+ij b/v+jv'] of JV elements 
of degree < k (because b^ G {x m+ i)k for JV' < i < JV + JV') such that 

} k = span(S) + (/ lr .,/ s ) w n%,... 

(here we use the notation of Assumption 13. 21 again). Thus any polynomial of degree 
< k can be rewritten, modulo (/i, . . . , f s )k+i, as a linear combination of elements 
in B of degree < k. As B contains 1 since b^, +1 = a^+i) this shows that B is a 
generating set of A. As .4. is of dimension JV, i? is in fact a basis of A, and thus 
S := k + 1 and fc satisfy the conditions in fl}. □ 

Remark 3.6. Note that in general 

(/i,...,/.> nK[x] d ^(h,...,f s ) d , 

where (/i, . . . , / s )d was defined in Inequality can happen when the system has 
a root at infinity, for example, if /i = x + 1, /2 = x i/ien (/i, /a) fl K[x]o = K 
&ui (/i,/2)o = {0}- However, using the homogenization fi,...,f^, the degree d 
part of the homogenized ideal is always equal to the space spanned by the multiples 
°f fii ■ ■ ■ > fs °f degree d. The above example also demonstrates that dim*4 is not 
always the same as dimK[x]d/(/i, ■ • ■ , f s )d even for large d, because above dim A = 
but dimK[x, y) d /{h,h)d = 1 for all d>0. 

Definition 3.7. Let JV := dimK(-4) and fix B = \b\, . . . , 6/v] a monomial basis for 
A such that dcg(6i) < k for all i = 1, . . . , JV. We define D to be the maximum 
degree of the monomials in B. Thus D < k < S. 
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Next we will define Sylvester and Macaulay type resultant matrices for /i, . . . f s . 
Definition 3.8. Define 

A :=max(<5- 1,2.0) 

where d and D are defined in Assumption \3. 2\ and Definition \3.7\ 
Let Syl A+1 (f) be the transpose matrix of the linear map 

(3) ®K[x]a- 4+ , — > K[x] A +i 

i 

s 

(gi,...,g„) i-> ^/ift 
i=i 

written in the monomial bases. So, in our notation, Syl A+1 (f) will have rows which 
correspond to all polynomials fix" of degree at most A. 

Let MacA (f ) be the matrix with rows corresponding to a basis of (fi , . . . , f s ) a+i H 
HC[x]aj obtained by eliminating coefficients of terms of degree A + 1 in the matrix 
Syl^ +1 (f) using Gaussian elimination, and then taking a maximal linearly inde- 
pendent set among the eliminated rows. 

Remark 3.9. In the case where s — m, for generic f with no roots at infinity, we 
can directly construct MacA(f) by taking the restriction of the map f3J) to 

m 

05,(A)^%] A 

i=l 

where <S,(A) = span{x Q : \a\ < A — dj, Vj < i, a,j < dj}. 

Here MacA(f) is a submatrix of the classical Macaulay matrix of the homoge- 
nization off and some fm+i> where f^+i * s an U homogeneous polynomial of degree 
A — 5: we only take the rows corresponding to the polynomials in f . Since the 
Macaulay matrix is generically non-singular, 
MacA(f) will also be generically full rank. 

Note that with our assumption that f\ , . . . , f m has no roots at infinity, we have 
that MacA(f) has column corank dirndl = YiT^i ^i- 

Since A > S — 1, by Assumption 13.21 and Theorem 13. 51 the corank of MacA(f) = 
N, where N is the dimension of A. Also, we can assume that the first columns of 
MacA(f) correspond to a basis B of A. 

Fix an element 

y=[y a :ae N"\ \a\ < A] T 
of the nullspace Null(Mac A (f)), i.e. Mac A (f) ■ y = 0. 

Definition 3.10. Let B = [b\, ...,6jv] be the basis of A as above, consisting of 
monomials of degree at most D. Using y we can define A y G by A y (g) := 
J2x a £B Va9a, where g = X) x Q es 5a x ° e -A. Note that every A ^ A* can be defined 
as A y for some y G Null(MacA(f)) or more generally with an element o/K[x]* 
which vanishes on the ideal X. 

Define the moment matrix 9Jls(y) to be the N x TV matrix given by 

9tts(y) = [y a +i3) a ,i3, 

where a and (3 run through the exponents of the monomials in B. Note that WIb is 
only a submatrix of the usual notion of moment matrix, see for example . 
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For p £ A, we define the linear function p ■ A £ A* as p ■ A(g) := A(pg) for all 
g£A. 

Remark 3.11. // one considers a linear function A on A, such that the bilinear 
form (x,y) i— > A(xy) is nondegenerate on A, then the moment matrix corresponding 
to this A will be the one whose (i,j)-th entry is just A(6,&j). Moreover, for g,h £ A 

A y {gh) = coeff B (<?) T • 97t B (y) • cocS B (h) 

where coe&sip) denotes the vector of coefficients of p £ A in the basis B. 

The following proposition is a simple corollary of [351 Prop 3.3 and Cor. 3.1]. 

Proposition 3.12. Let y be a random element of the vector space Null(MacA(f))- 
With high probability, 97ts(y) is non-singular. 

Remark 3.13. Using the above proposition, one can detect whether the algebra A 
is not Gorenstein with high probability by simply computing the rank of s HIb (y) for 
(perhaps several) random elements y in Null(MacA(f))- 

Remark 3.14. By [35J Theorem 2.6 and Lemma 3.2] one can extend y toy £ K Nm 
such that the infinite moment matrix 9Jl(y) := [y Q +/3]a,/36N m has the same rank as 
9Jte(y) and the columns o/9Jt(y) vanish on all the elements of the ideal X. 

Next we define a basis dual to B = [bi, . . . , 6jv] with respect to the moment 
matrix 9Jts(y). Using this dual basis we also define a polynomial J which is in some 
sense a generalization of the Jacobian of a well-constrained polynomial system. 

Definition 3.15. From now on we fix y £ Null (Mac a (f )) such that 2Ts(y) is 
invertible and we will denote by A the corresponding element A y £ A* . We define 

Let b* := Ylj=i c ji°j- Then [b*, . . . , b* N ] corresponds to the columns of the inverse 
matrix SJT^^y) and they also form a basis for A. Note that we have A(bib*) = 1, 
if i = j, and otherwise. 

Define the generalized Jacobian by 

N 

(4) J:=Y^bib$ modi 

i=l 

expressed in the basis B = [b±, . . . , 6/v] of A. 

Remark 3.16. Note that since y^.j bib* has degree at most 2D, and A > 2D, we 
can use MaCA(f) to find its reduced form, which is J. Because of this reduction, 
we have that deg( J) < D < 5. 

Note that the notion of generalized Jacobian was also introduced in [5] . Its name 
come from the fact that if s — m and if A is the so called residue (c.f. |17| ). then 
^2iLi °ib* = J is the Jacobian of /i, . . . , f m . 

We now recall the definition of the multiplication matrices and the matrix of 
traces as presented in [2"B"] . 
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Definition 3.17. Let p £ A. The multiplication matrix M p is the transpose of the 
matrix of the multiplication map 

Mp-.A — > A 

9 !->• P9 

written in the basis B. 

The matrix of traces is the N x N symmetric matrix: 

T=[Tr(b i b j )]% =1 

where Tr{pq) := Tr{M pq ), M pq is the multiplication matrix of pq as an element in 
A in terms of the basis B = [b±, . . . , 6at] and Tr indicates the trace of a matrix. 

The next results relate the multiplication by J matrix to the matrix of traces T. 

Proposition 3.18. Let Mj be the multiplication matrix of J with respect to the 
basis B . We then have that 

Mj = [Trihb*)}"^. 
Proof. Let A 6 A* be as in Definition 13.151 For any h £ A we have that 

N N 
N 

=> hb l = J2 A (h h b*)b,j M h [j, i] = A(h b, b*) 

N N 

=> Tr(h) = J2Mhbib*) = A(hY,bib*). 

j=l i=l 

Since J = Y$Li in A, we have Tr(h) = A(hJ). Therefore 
Mj\j ) i]=A(Jb i b*)=Tr(b i b*) 

□ 

Corollary 3.19. 

Mj-fm B (y) = [Tr(b i b j )}» j=1 =T, 
or equivalently J ■ A = Tr in A* . 

Proof. The coefficients of b* in the basis B = [6i, . . - , bjv"] are the columns of 
37l^ 1 (y), which implies that 

Mj = [Tribib*)}"^ = [Tr(&i&i)]& = i ■ Tl B \y). 

Therefore we have that Mj ■ M B {y) = [Tr(bibj)]F j=1 . □ 

Finally, we prove that the matrix of traces T can be computed directly from the 
Macaulay matrix of /i, . . . , f s and J, without using the multiplication matrix Mj. 
First we need a lemma. 

Lemma 3.20. There exists a unique matrix 5Hs(y) of size |Mon<(A) — B\ x \B\ 
such that 
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Mac A (f) 




= 



Proof. By our assumption that the first columns of MacA(f) correspond to B we 
have 




where the columns of B are indexed by the monomials in B. Note here that by 
Definition 13.81 and Assumption 13.21 the rows of MacA(f) span Ia+i n K[x]a, and 
the monomials in B span the factor space K[x]a/(2a+i n IK[x]a). These together 
imply that the (square) submatrix A is invertiblc. 
Then 



Id 



NxN 



A^B 



= 



which implies that 



Mac A (f) 



where *K B (y) = -A^B • 9Jt B (y). 




□ 



By construction, the column of OTs(y) indexed by bj 6 B corresponds to the 
values of bj ■ A € A* on £>i, . . . , bjy. The same column in EHs(y) corresponds to the 
values of bj ■ A on the complementary set of monomials of Mon<(A). The column 
in the stacked matrix corresponds to the value of bj ■ A on all the monomials in 
Mon<(A). To evaluate bj ■ A(p) for a polynomial p of degree < A, we simply 
compute the inner product of the coefficient vector of p with this column. 

Definition 3.21. Let B = [bi, . . . , 6jv] be the basis of A as above, and let P € K[x] 
be a polynomial of degree at most D. 

Define Syl B (P) to be the matrix with rows corresponding to the coefficients of the 
polynomials (b\P), . . . , (b^P) in the monomial basis Mon<(A) (we use here that 
deg(6i) < D, thus deg(6 4 P) < 2D < A). 

Furthermore, we assume that the monomials corresponding to the columns of 
Syl B (P) are in the same order as the monomials corresponding to the columns of 
Mac A (f). 
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Theorem 3.22. 




= [TrQbibj)]?^ 



Proof. Since the j-th column of the matrix 




represents the values of bj ■ A on all the monomials of degree less than or equal to 
A, and the i-th. row of Syl B ( J) is the coefficient vector of bi J, we have 




= [{b r A)(b i J)]? j=1 



[A(J^)fi = i 



□ 



We can now describe the algorithm to compute a set of multiplication matrices 
M Xi , i = 1, ... ,to of the radical \fl of T with respect to a basis of K[x]/v^- To 
prove that the algorithm below is correct we need the following result from |26l 
Proposition 8.3] which is the consequence of the fact that the kernel of the matrix 
of traces corresponds to the radical of A: 

Proposition 3.23. Let T be a maximal non-singular submatrix of the matrix of 
traces T. Let r be the rank of T, and B := [6^, . . . ,bi r ] be the monomials corre- 
sponding to the columns of T . Then B is a basis of the algebra K[x]/Vx and for 
each k = 1, . . . , m, the solution M x 



L Xk of the linear matrix equation 



TM Xk = T Xk 



is the multiplication matrix of for \/X with respect to B. Here T Xk is the r x r 
submatrix of [Tr(xkbibj)]f i j =1 with the same row and column indices as in T . 



Algorithm 3.24 (Radical ideal using Macaulay matrices and traces). 
Input; f = [fi, ■ ■ . , f s ] £ K[x] of degrees d±,...,d s generating an ideal X and 
5 > such that for k := S — 1 they satisfy the conditions in A ssumption \3.2i An 
optional input is D < S, which by default is set to be 8. 
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Output: A basis B for the factor algebra K[x]/vT and a set of multiplication 
matrices {M Xi \i = 1, . . . , to} of \/T with respect to the basis B. 



(1) Compute MacA(f) for A := max(2_D, 5 — 1) as in Definition \3.8[ 

(2) Compute a basis B of K[x]A/((f) a+i H K[x]a) such that the polynomials 
in B have degrees at most D. Let B = [b±, . . . , bff]. 

(3) Compute a random combination y of the elements of a basis of Null(Mac^(f)) . 

(4) Compute the moment matrix 971s (y) defined in Definition \3.10\ and ^s(y) 
defined in Lemma \3.20\ 

(5) Compute 9Jl^ 1 (y) and the basis [&*,..., b^] defined in DeHnition \3.15[ 



(6) Compute J = modT using MacA(f)- 

(7) Compute Syl s (J) and Syl B (xkJ) for k = 1, . . 

(8) Compute 



, to defined in Definition 



T 



T'ib.b,)^ 



and 



T Xk :=[Tr{x k hb 3 )]» 




Syl B {x k J) 



for k = 1, . . . , m. 



(9) Compute T, a maximal non-singular submatrix of T . Let r be the rank of 
T, and B := , . . . , 6,: r ] be the monomials corresponding to the columns 
off. 

(10) For each k = 1, . . . ,m solve the linear matrix equation TM Xk = T Xk , where 
T Xk is the submatrix ofT Xk with the same row and column indices as in T. 

Remark 3.25. Since the bound given in Theorem Iff. 51 might be too high, it seems 
reasonable to design the algorithm in an iterative fashion, similarly to the algo- 
rithms in [35., 36, 46], in order to avoid nullspace computations for large matrices. 
The bottleneck of our algorithm is doing computations with MacA(f), since its size 
exponentially increases as A increases. 

Remark 3.26. Note that if s = to then we can use the conventional Jacobian of 
fx, . . . , f m in the place of J, and any |Mon<(A)| x \B\ matrix X such that it has 
full rank and MacA(f) ■ X = in the place of 
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Even though this way we will not get matrices of traces, a system of multiplication 
matrices of the radical \fl can still be recovered: if Q denotes a maximal non- 
singular submatrix o/Syl B (J) • X, and Q Xk is the submatrix of Syl B (xk J) • X with 
the same row and column indices as in Q, then the solution M Xk of the linear matrix 
equation QM Xk = Q Xk gives the same multiplication matrix of \/X w.r.t. the same 
basis B as the above Algorithm. 

Remark 3.27. As M Xk is the transpose matrix of multiplication by Xk modulo the 
radical ideal \fX, its eigenvectors are (up to a non-zero scalar) the evaluation at 
the roots Q of X (see [381 117] for more details). The vector which represents this 
evaluation at Q in the dual space A* is the vector of values of [b±, . . . , bj^] at £. To 
obtain these vectors, we solve the generalized eigenvalue problem (T Xk — zT) w = 
and compute v = Tw. The vectors v will be of the form [&i(C), ■ • ■ , °n(C)) f or C a 
root of T. If b\ — 1, 62 = xi, . . . , bm+i = x m , we can read directly the coordinates 
of £ from this vector. 

3.2. The Non-Gorenstein Case. We will now consider the case where A is not 
Gorenstein. The main idea of the algorithm is the same as in the Gorenstein case, 
except we will obtain as an output a matrix of traces with respect to an algebra B 
which is a maximal Gorenstein factor of A. This will still allow us to compute the 
multiplication matrices of the radical of T since the maximal non-singular submatrix 
of the trace matrix corresponding to B is the same as that of the trace matrix of 
A. First we will need some results to define a maximal Gorenstein factor B of A 
from a random element of the nullspace of MacA(f). 

Let K be an arbitrary algebraically closed field. All algebras we consider will 
be finite dimensional commutative K-algcbras. A local K-algebra here is an K- 
algebra B, with unique maximal ideal (which we denote by M.) such that Bj M. is 
isomorphic to K. Note that due to the fact that IK is algebraically closed, no other 
residue class field is possible. 

Definition 3.28. Fix y <G Null(MacA(f)) and let A := A y e A* defined as in 
Definition \3.1(A We define 

11(A) := {a e A, A(ab) = holds for all b 6 A}. 

Note that 1Z(A) = iff A(xy) is a nondegenerate bilinear form on A. Also, an 
easy calculation shows that 1Z(A) is an ideal in A. 
Define 

B :=A/K(A). 

First we need the following technical lemmas. 

Lemma 3.29. 1Z(A) is an ideal of A, in fact, it is the largest ideal of A which is 
in ker(A). 

Here ker(A) is the set of elements a G A, such that A(a) = 0. 

Proof. An easy calculation shows that H(A) is an ideal. Clearly it is in ker(A). 
Conversely, if J is an ideal of A which is in ker(A), then J is in 7?.(A). Note that 
the sum of ideals is an ideal again, hence there exists a unique largest ideal of A 
which is a subset of ker(A). □ 

Lemma 3.30. (Structure Theorem on Artininan Rings, specialized to our setting, 
see Atiyah-MacDonald) A 'K-algebra B is the direct product of finitely many local 
WL-algebras. 
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Lemma 3.31. Suppose that the "K-algebra B is the direct product of the local K- 
algebras Bi, (i = 1, . . . , k). Then B is Gorenstein iff all the Bi are Gorenstein. 

Proof. Note that the Bi can be viewed as ideals of B, moreover B is the direct sum 
of these (as K-subspaces) . 

Let A be a linear function on B such that the form A(xy) is nondegenerate. Then 
the resctriction A.; of A on Bi will define a nondegenerate bilinear form on Bi. Let 
a be a nonzero element of Bi. Then there exists a 6 £ B such that A(ab) ^ 0. Write 
now b as b = ^ bj , with bj £ Bj for ( j = 1 , . . . , k) . Note that we have abj = if 
j 7^ i, hence 

0^ A(ab) = A(ah) = A t (ab % ), 

hence Ai(xy) is nondegenerate. 

Conversely, assume that we have linear functions A^ : Bi — > K such that the 
form Ai(xy) is nondegenerate on Bi. We define A as follows. Let a £ B with 
a = ai + ■ • • + dfc, where ai € Bi. Note that a uniquely determines the a, , and the 
map a i— > otj is an K-algebra morphism from B to Bi. This implies that A(a) := 
Ai(ai) + - • -+Afc(afe) is a correct definition and A is a linear function on B. Moreover, 
it is easily seen that A(xy) is nondegenerate. □ 

Lemma 3.32. A induces a linear function A' on the factor B = A/1Z(A). For the 
function A' on B we have that IZ(A') = ; hence B is Gorenstein. 

Proof. We set A' (a + 71(A)) = A(a). It is routine to check that A' is a correctly 
defined linear function on B. Suppose that a + TZ(A) is in TZ(A'). Then A(ac) = 
holds for every c £ A, hence a £ 71(A) and therefore a + 71(A) = in B. □ 

Lemma 3.33. Every Gorenstein factor B of an "K-algebra A can be obtained via a 
linear function A on A, as outlined by Lemma \3. 32i 

Proof. Let B = A/ J be a Gorenstein factor of A. Let A' be a linear function 
on B with TZ(A') = 0. We define A on A as follows. Let A(a) := A' (a + J) for 
a £ A. Clearly A will be a linear function on A. Let a £ A such that A(ab) = 
for every b £ A. Then = A(ab) = A'(ab + J) = A'((a + J)(b + J)), giving that 
a + J £ Ti(A'), hence a £ J . This shows that 7?.(A) C J . The reverse containment 
is immediate, therefore 7Z(A) = J . Now one can directly check that A' is obtained 
from A via the factor construction of Lemma 13.321 □ 

Lemma 3.34. Let A be an ^K-algebra. Then the maximal dimensional Gorenstein 
factors of A are obtained as follows: from every direct factor Ai (this is local) factor 
out an ideal Ji so that Ai j Ji is a maximal Gorenstein factor of Ai . 

Proof. Any ideal J of A is the direct sum of the ideals Jj = Aj fl J . On the other 
hand Aj J is Gorenstein iff the local factors Ai/ ' Ji are (Lemma 13.311) . □ 

The following Theorem shows that we can get maximal Gorenstein factors of A 
from random linear forms on A with high probability, similarly as in the Gorenstein 
case. 

Theorem 3.35. The maximal Gorenstein factors of A can be obtained with high 
probability as B := A/1Z(A), where A is a random linear function on A. 
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Proof. Let a\, . . . , a m be a basis of A over K. Then for a linear function A : A — ► IK 
the dimension of A/1Z(A) is the rank of the matrix [A(ojOj)]. Thus, maximal 
Gorenstein factors are obtained if the rank of the matrix is maximal (as A ranges 
over the K-dual of A). If one fixes a dual basis of A, and writes A as a linear 
combination of these basis functions, then, the entries of the matrix [A(aiOj)] will 
be linear polynomials of the coordinates 71 , . . . , j m of A. Now consider a A which 
achieves the maximal rank k, and consider a corresponding k x k minor of the 
matrix that has a nonzero determinant. This determinant is not identically zero, 
as a function of the 7^, hence it will be nonzero on a Zariski open set. The linear 
functions corresponding to the points in this set will define maximal Gorenstein 
factors. □ 

We now show that any maximal Gorenstein factor will allow us to compute the 
radical of A. 

Theorem 3.36. Assume that A is such that the corresponding bilinear form on A 
has maximal rank. Then 

K{A) C Rad(^). 

Proof. By Lemma T3.341 we can assume that A is local. By Lemma f3. 291 TZ(A) is an 
ideal. Since ^4/Rad(^4) = K is Gorenstein and B is the maximal Gorenstein factor 
of A we have that 7£(A) 7^ A. Therefore 72.(A) is a subset of the unique maximal 
ideal M = Rad(A). □ 

Using the previous results, we are now ready to define the main ingredients of 
our algorithm in the non-Gorenstein case, which are analogous to the Gorenstein 
case except that instead of working in A we are going to work in B. We will use 
the notation of the previous subsection. 

We can obtain a basis for B as follows. Let B = [bi, ... , 6jv] be a basis for A 
and y 6 Null(MacA(f)) such that the moment matrix SCls(y) has maximal rank. 
Since the columns of OTs(y) correspond to B, taking a maximal nonsingular minor 
of this matrix will define a subset B a = [b ai , . . . , b Ur ] of B corresponding to the 
columns of this submatrix. Then B a will form a basis for B as we prove in the 
following proposition. 

Proposition 3.37. B a = [b ai , . . . , b ar ] forms a basis for B. 

Proof. Consider the moment matrix 9Jl_e(y) with columns corresponding to B and 
let r be its rank. Since B a corresponds to a set of basic columns of 97tg(y), there 
exists a basis Vx, ... , vjv_ r for Null(9Jts(y)) which can be extended to a basis of 
K N by adding the unit vectors e ai := [S ai ,j}jLi for i — 1, . . . , r. 

Let Vi be the element of A obtained by taking the linear combination of 61 , . . . , 6jv 
corresponding to the coordinates of Vj for i = 1, . . . , N — r. Then it is easy to see 
that Vi € 1Z{A). Thus the elements of B a correspond to a basis of A/H(A) = B. □ 

Here we need to define the moment matrix. 

Definition 3.38. Let B a be defined as above. Define 

k 
1=1 

similarly as 
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9ttB„(y) 



9Wy) 



Using the following theorem we get that the maximal nonsingular minor of the 
smaller trace matrix [Tr(b ai b aj )]*j =1 suffices to compute the radical. 

Theorem 3.40. Let B — [b±, . . . , 6jv] be a basis for A and B a = [b ai , . . . , b ar ] a 
basis for B, where ol\ < ■ ■ ■ < a r are in {1, ... , N}. Then 

rank[Tr(6 Qi 6 Qj )][ J=1 = rank[Tr(6 i 6 i )]^ =1 . 

(As before, by a slight abuse of notation, we use the same notation for elements in 
A, in B, and their common preimages in K[x].j 

Proof. This follows from Theorem 13.361 and the fact that the rank of the trace 
matrix is dim.A/Rad(.A) = dimS/Rad(S). □ 

4. Affine Complete Intersection Ideals 

4.1. Univariate Case. In this section we will follow the work of Mourrain and 
Pan [40]. We start by defining the univariate Bezout matrix. 

Definition 4.1. Let f,g G K[x] be two univariate polynomials such that degg < 
deg f = d, and let y be a new variable. Then the Bezoutian O f t3 of f and g is the 
polynomial 

" 0<i,j<d-l 

The Bezout Matrix Bf g of f and g is the d x d matrix 

[Bf,g]ij — Cij ■ 

We will need the following definition: 

Definition 4.2. The Horner basis for the polynomial f is the set {Hd-i, ■ ■ ■ , Hq} 
with 

Hi(x) = a i+ \ + h a d x d ~ l ~ 1 for i = 0, . . . , d — 1. 

Note that in terms of the Horner polynomials, we have that 

d-l 

<B f , 1 (x,y) = J2x i H i (y). 

The following theorem connects the Bezoutian of / and its derivative /' with 
the matrix of traces of / with respect to the Horner basis. 

Theorem 4.3. 

B fJ , = [Tr(HiHj)} d J* , 
where [Tr(HiHj)] d ~]^ is the matrix of traces of f in the Horner basis (see eg. [7], 



Theorem 3.39. [Tr(b ai b aj )] 



k 
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Theorem 14.31 implies that using Dickson's Lemma one can compute the square- 
free factor of / by simply computing the kernel of Bfji. It's natural to ask how 
our method based on Dickson's lemma relates to computing the square-free fac- 
tor of / via computing gcd ^ ■ The following proposition shows that computing 
//gcd(/, /') to get the square- free factor using the Bezout matrix is computation- 
ally equivalent to using Dickson's Lemma. 

Proposition 4.4. The smallest degree polynomial of the form 2i=o r iHi{x) such 
that [ro, . . . ,r^_i] T is in the kernel of Bfji, is equal to //gcd(/, /'). 

4.2. Multivariate Case. For the multivariate case we will first define the mul- 
tivariate analogue of the Bezout matrix (also referred to as Dixon matrix in the 
literature). The papers [5] and [57] are good references for the Bezout (Dixon) 
matrix described below. 

Definition 4.5. Let 

f := [/i,.,/ m ]6K[ I i,.,x m ] m 

and consider an additional polynomial /o € K[xi, . . . , x m \. We use the notation 
x = [xi, . . . , x m ], y = [yi, . . . , y m ] and 

X = [xi, . . . ,x m ],Xi = [yi,x 2 , ■ ■ -,x m ], ■ ■ ■ ,X m = [yi, y 2 , ■ ■ ■ , y m }- 

The Bezoutian of the system [/o, fx , . . . , f n ] , denoted by 53 f , is a polynomial in the 
variables x and y defined as follows: 

fo(Xo) 

/m(*l) 



930§Xx,y) :=dct 



/o(*o) 


- /o(*i) 


fo(X m -i 


-fo(Xm) 


Xi 

h(Xo) 


- h(Xi) 


Xrn 

fl(X m -i 


- y m 

-fl(Xm) 


Xl 


- yi 


Xm 


~ Vm 


fm{Xo) 


— fm(Xi) 


fm- (X m — 1 


~ fm(X m ) 


Xl 


- yi 


Xm 


- y m 



The Bezout matrix of the system [/o, /i, • ■ ■ , fm], denoted by Bf , is the coefficient 
matrix of the Bezoutian, i. e. if we write 

<8 /0 (x,y)= c Qi/3 (/o)xV 

where E and E' are subsets of W n and c a ,p(fo) € K, then the Bezout matrix of 
[fo, fi ■ • ■ , fm] is the \E'\ x \E\ matrix 

We denote by 93 * the map 

(6) <B% : (\p) 0eE , ^ J2 CoAfo)**^ 

aEE, 0EE' 

and by Q3j o the map 

(7) : (A a ) agB -> ^Afo)X a y f3 - 

Our goal is to compute a matrix of traces for the system f = [f±, . . . , f m ] from 
the Bezout matrix Bj of the system [ J, f\ , . . . , f m ] analogously to the univariate 
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case, where J is the Jacobian of /i, . . . , f m - As we mentioned in the introduction, 
in general the Bezout matrix Bj is not a matrix of traces, which can be easily seen 
by comparing sizes. However, to obtain a matrix of traces of f one can define a 
reduced version of the Bezoutian and the Bezout matrix as follows. 

Definition 4.6. Let f = [/i, . . . , f rn ] G K[xi, . . . , x m ] m and assume that the factor 
algebra *4(x) := K[x]/X has dimension TV over K, where I is the ideal generated 
by the polynomials in f. For some /q G K[x] let Q5/ (x, y) be the Bezoutian of the 
system /o, /i, . . ■ , f m . Let B = [b±, . . . , bpj] and B' = [b[, . . . b' N ] be bases for A(x) 
and A(y), respectively. Then we can uniquely write 

® /o (x,y) = E Pb,v(fo)b(x)b'(y) + F( X ,y) 

b£B,b'£B' 

where F(x, y) G (I(x),/(y)) and Pb,b'(fo) G IK. VFe define the reduced Bezoutian 
25 / wzi/i respect to the bases B and B' as 

»/b~ E /?M'(/o)6(x)6'(y) 

&G-B,b'eB' 

and i/ie reduced Bezout matrix Bf with respect to the bases B and B' to be the 
TV x TV matrix 

B h B : = (/ 3 b,b'(/o)) 6e ij.6'es' • 
We are going to use the following theorem [T^] : 

Theorem 4.7. There exists (dual) bases 9 := [#i(x), . . . , 6*jv(x)] and il := [u>i(y), . . . , wjv(y)] 
of A(x) and A(y) such that for all polynomial f G K[x], we have 

»/ : = E <M/) Wj (y) + F(x, y), 
wii/i -F(x, y) G /(x) <g) /(y) and such that 

i 

In [40] the following expression was given for the reduced Bezout matrix Bj of 
the system [J, f±, . . . , f m ] in terms of a matrix of traces of f±, . . . , f m : 

Theorem 4.8. Let f, A(x), A(y) and the bases = (6*;), f2 = (cjj) 6e as m 
Definition ^. 6\ and Theore m\4.7\ Let J be the Jacobian of fi, . . . , f m , and consider 
the reduced Bezout matrix 23 j with respect to the bases and 0. Then 

B*j' B = [Tr{u i u j )]» j=1 . 

Using the relation B/ = /3j'^/3f' Q , we deduce that 

B?' n = [Tr( Wi ^ =lJ 

so that A = [Ai] G kcrSj'° iff 

JV 

Tr(w i (5^A i e i )) = 0, i = l,..., TV, 
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or equivalently iff 

N 

r:=X>0 3 (x)eVI. 

i=i 

Because of the block diagonal form of the Bczoutian matrices in a common basis 
(Theorem 14. 7p . we deduce that if A is an element of ker(Q5j) then 

Q5?(A) = r(x)+/i(x), 

where h £ /, Q5i is the Bezoutian matrix of 1 in the (monomial) bases (x q )q, 6 £, 
{yp)peE', ® i is the corresponding map defined in ([5]), and im(25^) is the space 
generated by the coefficient vectors with respect to (y^)p^E' of the polynomials in 
the image of the map 25^ (see 0). Then, we have the following theorem: 

Theorem 4.9. Using the previous notation we have that 

v7=B?(ker(»5))+J(x). 

Proof. Because of the block diagonal form of the Bezoutian matrices in a common 
basis (Theorem 14. 7 j) and the previous discussion, we deduce that if A is an element 
of ker(25j) then 

<8 X 7 (A) ^A^^A^fy))!^) = 

where F A (x) £ I. By the previous discussions, 

^■(x)A(^-(y))=!B x (A)eV7. 

j 

□ 

Note that the role of x and y can be exchanged in this theorem. Note also that 
kcr(Q5j) can be replaced by ker(Q5j) n ker(Q5^) ± in this theorem. 

A question that remains is how to compute the multiplication matrices M X1 , . . . , M. 
of the radical \fl. 

In order to compute the reduced Bczout matrix B / of the system [f\ , . . . , f m ] 
with respect to some bases [oi, . . . , a r ] of A(x) and . . . , b r ] of A(y), it is sufficient 
to find expressions of the form 

x a = ^c 6 fe(x) + F(x), for all a £ E and 
y 13 = Y, c b>b'{y) + G(y), for all (3 £ E' 

b'eB' 

where F,G £ I, Cb, cy £ IK and E and E' were defined in Definition 14.51 Define 

V := (x a |a £ E) and W := (y p \f3 £ E'). 

Assuming that 6j £ V and b' { £ W for i = 1, . . . , I, the task is to find enough linear 
combinations of the monomials corresponding to the rows and the columns of the 
Bczout matrices which belong to the ideal I. 

We follow the approach described in [39] , where it was shown that the computa- 
tion of the Bezoutians %$ Xi of the system [xi, fi, . . . , f m ], for i = 1, . . . , m, as well 
as the Bezoutian 05 1 of the system [1, /i, . . . , f m ], gives sufficient information of the 
structure of / in order to find the reduced Bezout matrix Bf for any f £ K[x]. In 
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order to get the structure of yfl we simply have to add the polynomials in x (resp. 
y), obtained from Theorem l4.9l 

Here we describe a summary of this method. First notice that 

x^&x - ©x, € J(x) and la®! - B w G /(y) i = 1, . . . , m. 

The initial step of the method is to obtain ideal elements in \/T(x) (resp. 
\/T(x)) which are in V (resp. W) from a:j5Bi — Q3 Xj (resp. y^QSi — 58jk) and also 
from 03* (ker(*Bj)) ( resp. <B( (ker(Sj))). The elements in V7(x)n V obtained by 
the initial step are denote Kq, and the ones in y/l(y) f~l W are denoted by Hq. 

For any vector space K C R, we denote by K + , the vector space K + = K + 
X\K + - ■ - + x m K. The notation means n iterations of the operator +, starting 
from K . 

To prove that we get the quotient structure by the radical ideal \/T, we will 
assume that V is connected to 1, that is. V contains 1 and for any v G V — (1), there 
exists I > such that v G span(l)W and v = Uo+X^L i x i v i with G span(l)" -1 ]nV 
for i = 0, . . . , m. 

In order to obtain additional ideal elements, the following steps are used [55] : 
Saturation step: Finds new ideal elements by multiplying the already com- 
puted ideal elements by the variables Xi for all i = 1, . . . , m. 
Column reduction step: Finds new bases for the vector spaces V and W 
such that the new basis for V contains previously computed elements in 
VT(x) n V, and also that the Bezout matrix B\, written in terms of these 
new bases, has a lower block triangular structure. By writing the matrices 
B Xi in terms of the new bases for V and W , one can obtain new elements 
in Vl(x) n V. 

Diagonalization step: After the column reduction step one can transform 
Bi into a block diagonal form which, by repeating the same transformation 
on the matrices B Xi , can possibly reveal new ideal elements. 

Row reduction step: Same as the column reduction step, with the roles of 
x and y interchanged. 
They are used in the following iterative algorithm: 

Algorithm 4.10 (Radical of an affine complete intersection). 
INPUT: f = . . . , f m ] G K[x] generating an ideal I, which have a finit number 
of complex roots. 

Output: M Xl , . . . , M x a system of multiplication matrices for the radical ideal 
sfl. 

• Compute the Bezoutian matrices B\ , B Xl , . . . , B Xm of 1 , x\ , . . . , x m and /i , . . . , /, 
and Bj. 

• Using the initial step, define K := Kq; H := Hq; notsat := true. 

• While notsat 

— Apply the saturation step on K and H ; 

— Apply the column reduction step; 

— Apply the diagonalisation step; 

— Apply the row reduction step; 

— If this extends strictly K ; or H , then 

let notsat := true, otherwise let notsat := false. 
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• Return 

(8) M Xi :=N^N Xi , i = l,...,m 

where N Xi is the matrix reduced from B Xi , at the end of the loop. 

The loop terminantes because the size of the matrices are decreasing. If it ends 
with matrices of non-zero size, iVx is necessarily invcrtiblc. At the initial step and 
all along the computation, we have K c \fl(x),H C VT(y), since the different 
steps are valid modulo / c vT. We denote by [ai, . . . , a r ] (resp. [bi, . . . ,b r ]) the 
linearly independent polynomials indexing the rows (resp. columns) of Mi and A 
(resp. B) the vector space they span. By construction, the vector space V (resp. 
W) decomposes as V = A + span(X) (resp. W = B + span(iJ)) where K and H 
are the sets of relations in \fl updated in the reduction steps during the algorithm. 
We complete a,i,i = 1, . . . , r (resp. bi, i = 1, . . . , r) in a basis ax, ... , cl\e\ °f ^ (resp. 
bi, . . . , bifiii of W) with ai 6 K (resp. &j e i/) for i > r. 

We have the following theorem, which allow us to compute the radical of an 
affinc complete inetrsection, based on simple algebra tools: 

Theorem 4.11. Let fx,.-.,f m be as above. Upon termination, Algorithm \4-lC\ 
computes new bases [a\, . . . ,cl\e\] for V and [b±, . . . ,b\E'\] f or W, such that 

• [a\, . . . ,a r ] (resp. [b\, . . . ,b r ]) is a basis o/K[x]/v / /(x) (resp. K[y]/v / 7(y)J, 

• the ouput matrices Mi are the (resp. transpose) matrix of multiplication of 
Xi with respect to the basis [ax, . . . , a r ] (resp. [6i, . . . , b r \) for i = 1, . . .m, 

• a r+ i, . . . ,a\ E \ G \/T(x) n V and b r+x , . . . , b\ E <\ G VI(y) n W. 

Proof. The proofs of [33] [Lemma 5.3, Proposition 5.3-5.6] apply also here in order to 
show that the matrices M Xi commute. Moreover, by construction, K C V n \/T(x) 
such that V = A © span(iC), and K satisfies the following relations: fa £ (K) for 
i = 1, . . . , m (same proof as in [39] [Proposition 5.7]), (K) C VI by definition of 
the reduction steps, Bf(keiBJ) C (K ) C (K). Therefore, by Theorem [49] this 
shows that (K) = \fl. As in [39] [Theorem 5.8], we can assume that 1 £ A and 
that we have an exact sequence: 

-> (K) -> K[x] -> A -> 
/ ^ /(M)(l) 

where /(M) := f(M Xl , . . . , M lm ) is the linear operator of A, obtain by replacing 
the variable Xi by M Xi . As (K) = s/l, this shows that A ~ K[x]/-\// and the basis 
of A is a basis in K[x]/V7. □ 

5. Conclusion 

In an earlier work we gave an algorithm to compute matrices of traces and the 
radical of an ideal X which has finitely many projective common roots, none of 
them at infinity and its factor algebra is Gorenstein. The present paper considers 
an extension of the above algorithm which also works in the non-Gorenstein case 
and for systems which have roots at infinity, as well as an alternative method using 
Bezout matrices for the affine complete intersection case to compute the radical 

VT. 
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